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Abstract 

We have examined the frequencies of trapped two-armed (m = 2) nearly vertical 
oscillations (vertical p-mode) in vertically isothermal disks with toroidal magnetic 
fields. The magnetic fields are stratified so that the Alfven speed, ca, is constant in 
the vertical direction. We have particularly focused our attention on how frequencies 
of the fundamental mode {uj. = 0) and first overtone (rZr = 1) in the radial direction 
change with correlation, when the ratio c\/ changes, Cg being the isothermal acoustic 
speed. The results show that in the case where the oscillations are fundamental 
mode (ra = 1) in the vertical direction, the correlated frequency changes of the above- 
mentioned oscillations seem to well describe, with standard values of the mass and spin 
of the central sources, the frequency correlation of kHz QPOs observed in neutron-star 
X-ray binaries. 
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1. Introduction 

Many authors now think that the kilo-hertz quasi-periodic oscillations (kHz QPOs) ob- 
served in low-mass X-ray binaries (LMXBs) are some kinds of disk oscillations in a strong gravi- 
tational field. Understanding their origin is of importance, since it will become a promising tool 
to evaluate the mass and spin of the central sources as well as the disk structure surrounding the 
sources. Many models of kHz QPOs related to disk oscillations and waves have been proposed. 
None of them, however, can sufficiently well describe the observed frequency correlation of the 
twin kHz QPOs (Lin et al. 2010). 

As one of possible models of kHz QPOs, for example, we proposed a disk-oscillation 
model in deformed disks (Kato and Fukue 2006), based on the fact that in deformed disks a 
set of oscillations with positive and negative wave energies are excited resonantly (Kato 2004, 



2008a,b; Ferreira and Ogilvie 2008; Oktariani et al. 2010; Kato et al. 2011; Kato 2011b). This 
model qualitatively describe the frequency correlation of observed twin QPOs. However, there 
are some problems. The major problems are that the mass required to describe the observed 
frequency correlation of the twin QPOs of neutron-star X-ray binaries is as large as 2AMq and 
further that the calculated correlation deviates gradually from the observed one as frequency 
increases (Kato 2007; Lin et al. 2010). 

Other disk oscillation models we have are those in which trapped oscillations are directly 
connected to QPOs. In geometrically thin disks, there are many kinds of disk oscillation modes. 
In terms of the node number in the vertical direction, n, and the frequency in the corotating 
frame, these oscillation modes are classified into four, i.e., p-, g-, c-, and vertical p-modes (for 
details of the classification, see Kato 2001, Kato et al. 2008). Some oscillations of these modes 
are known to be trapped in the innermost region of relativistic disks, and thus they become a 
set of oscillations with discrete frequencies. These trapped oscillations have been studied by 
many authors, e.g., by Kato and Fukue (1980), Ortega-Rodrigues et al. (2002), Lai and Tsang 
(2009) for p-mode oscillations; by Okazaki et al. (1987), Nowak and Wagoner (1992), Perez 
et al. (1997), Fu and Lai (2009) for g-mode oscillastions; by Kato (1990), Silbergleit et al. 
(2001) for c-mode oscillations; and by Kato (2010, 2011a) for the vertical p-mode oscillations 
(see Wagoner 1999, Kato 2001, and Kato et al. 2008 for reviews of trapped oscillations). There 
seems, however, to be no quantitative studies examining whether a set of these oscillations 
can really describe the correlated frequency-changes of twin kHz QPOs. Here, we discuss this 
problem, by especially focusing on the vertical p-mode oscillations, using the results obtained 
by Kato (2011a). 

2. Trapped Vertical p-mode Oscillations in Magnetized Disks 

Here, we focus our attention on oscillations which have at least one node, n> 1, in the 
There are two oscillation modes with a given ra, when n > 1. One is the 
g-mode, the square of whose frequencies in a corotating frame, {u — mfi)^, is smaller than the 
square of the epicyclic frequency, k^, i.e., {u — mQy < k^, where u is the frequency of oscillations 
in the inertial frame, m is the azimuthal wavenumber and Q is the angular velocity of the disk 
rotation. The other mode of oscillations with > 1 is the vertical p-mode, the square of whose 
frequencies in the corotating frame is high in the sense that {u — mfi)^ > Q^, where Q± is the 
vertical epicyclic frequency and always larger than k. In this paper we consider the vertical 
p-mode oscillations of n = 1 and 2 with m = 2. 

To understand characteristics of trapping of the vertical p-mode oscillations, let us first 

^ When we mention the node number in the vertical direction, it is the number of node(s) of the density 
perturbation in the vertical direction. In the case of n = 1, the density perturbation has a node on the 
equatorial plane. It is noted that the node number associated with the vertical component of velocity is 
smaller than that of the density perturbation by one. 
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consider the case where the disk is vertically isothermal and has no global magnetic fields. 
Then, the propagation region of the vertical p-mode oscillations with u, m, and n is specified 
by (w — mf2)^ >nQ\ (e.g., see Kato 2001). This means that the propagation regions are divided 
into two radial regions specified by w < mQ — n^^'^^lj_ and uj > mf2 + n^/^r2_L§| The former shows 
that the propagation region of prograde oscillations, u > 0, with m = 2 and n < 4 is limited 
only in the innermost region of disks, and these oscillations are trapped there. That is, we can 
expect a discrete set of trapped oscillations. 

Now we consider the case where global toroidal magnetic fields exist in disks. Since the 
vertical p-mode oscillations are oscillations of nearly vertical motion, they are the fast mode of 
three MHD waves, when toroidal magnetic fields exist. Hence, the square of the frequency in 
the corotating frame becomes larger than n^l\ by the effects of magnetic restoring force, say 
{u — ■mQ±)'^ > {n + a)Q\, where a is a dimensionless quantity of the order of unity, depending 
on n and the strength of magnetic fields (see Kato 2011a for details). Hence, the propagation 
region near the inner region of the disks is changed to a; < mf2 — {n + a)^/^f2_|_. Since the right- 
hand side of this inequality is smaller than mf2 — n^^^fix, we expect that the eigen-frequency 
of the trapped oscillations become smaller than that in the case of no magnetic fields, if other 
parameters are fixed (see Kato 2011a). The frequency of the trapped oscillations depends also 
on the node number of oscillations in the radial direction (say, n,-), and decreases with increase 
of rij- if other parameters are fixed (see Kato 2011a). 

In this context, we examine in this paper whether the following picture can describe 
the observed frequency- correlation of twin kHz QPOs with reasonable parameter values. That 
is, we regard the fundamental mode of the vertical p-mode oscillations with no node (rij. = 0) 
in the radial direction as the upper kHz QPO, and the first overtone {riy- = 1) in the radial 
direction as the lower kHz QPO. Then, we examine whether the correlated time variation of 
the above two oscillations (tTt = and n,- = 1) due to time variation of magnetic fields can 
describe the correlated time variation of observed twin QPOs. It is noted that time variation 
of global magnetic fields will be natural, since they are amplified by winding due to differential 
rotation and decreases by dissipation due to reconnection. 

3. Calculation of Frequencies of Trapped Vertical p-mode Oscillations 

Frequencies and their parameter dependences of trapped vertical p-mode oscillations 
in disks with toroidal magnetic fields have been numerically calculated by Kato (2011a) by 
a perturbation method. Since we use here the results, we briefiy summarize the procedures 
adopted by Kato (2011a). 

^ Among oscillations of {ui — mfl)^ > 51^, one-armed (m — 1) one with one node {n — 1) in the vertical direction 
has a particular position. The oscillation is called c-mode (corrugation mode). It is nearly incompressible 
motions changing disk plane up and down with a corrugation pattern. This mode is outside of our interest 
in this paper. 
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The Newtonian formulation is adopted except that the effects of general relativity are 
taken into account when we condiser the radial distributions of fi(r), K,{r), and Q±{r). The 
angular velocity of disk rotation, Q{r), is approximated by the relativistic Keplerian one, since 
geometrically thin disks are considered. We adopt the cylindrical coordinates (r, cp, z), where 
the 2-axis is perpendicular to the unperturbed disk plane and the origin is at the disk center. 

3.1. Vertical Structure of Unperturbed Disks with Toroidal Magnetic Fields 

The unperturbed disks are assumed to be axisymmetric with toroidal magnetic fields. 
The fileds are assumed to be purely toroidal with no poloidal component: 

B,{r,z) = [0,Bo{r,z),0]. (1) 

We further assume that the gas is isothermal in the vertical direction and BQ{r,z) is distributed 
in such a way that the Alfven speed, ca, is constant in the vertical direction, i.e., (-Bg/47rpo)^^^ = 
const, in the vertical direction, where pQ{r,z) is density in the unperturbed disks. 

Then, the hydrostatic equilibrium in the vertical direction gives that pQ{r,z) and Bo{r,z) 
are distributed in the vertical direction as 

po(r,2;) =poo exp(^- and Boir.z) = Bqo exp(^- (2) 

where the scale heigth, H, is related to Cg, ca, and Q± by 

H\r) = (3) 

3.2. Calculations of Frequencies of Trapped Oscillations 

Small-amplitude perturbations are superposed on the equilibrium disks described above. 
The perturbations are taken to be proportional to exp[i{ut ~ rmp)]. Then, equations describing 
the perturbations are partial differential equations with respect to r and z [see equation (20) 
by Kato 2011a]. The equation is solved by being approximately decomposed into an equation 
describing the behaviour in the vertical direction and that in the radial direction. When doing 
so, we remember that we are now interested in the vertical p-mode oscillations, and thus among 
three velocity components of Ur, u^, and Uz, associated with the perturbations the main one is 
Uz- Furthermore, in the lowest order of approximations, the 2;-dependence of Uz can be expressed 
in terms of the Hermite polynomials Hn-iiz / H) , where n is a positive integer representing node 
number of oscillations in the vertical direction (Okazaki et al. 1987). Considering these facts, 
and using a perturbation method, we derive an ordinary differential equation with respect to r, 
which describes the behavior of perturbations in the radial direction. The equation is solved by 
a WKB method with relevant boundary conditions in the radial direction to obtain the eigen- 
frequency, u, of trapped oscillations. As an inner boundary condition to solve the equation 
we adopt that at the inner edge of the disk (which is taken at the radius of marginary stable 
circular orbit) Uz vanishes, and as the outer boundary condition we assume that the outside 
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of the trapped region (the evanescent region), Uz decreases outwards. The detailed procedures 
and the results on frequencies of trapped oscillations are given by Kato (2011a). 

The frequencies of trapped oscillation depend on various parameters specifying the disk 
structure (as well as the central sources) and the modes of oscillations. As the former types of 
parameters we have c^/Cg, M (mass of the central source), and a* (spin parameter). As the 
latter ones, we have m, n, and tt^. As mentioned before, we consider two-armed oscillations, 
i.e., 771 = 2. Concerning n and rij-, we restrict our attention in this paper to n = l and 2 (i.e., the 
fundamental and first overtone in the vertical direction), and nr = and 1 (i.e., the fundamental 
and first overtone in the radial direction). 

4. Frequency Changes and Their Correlation 

By using the procedures mentioned in the above section, Kato (2011a) numerically 
calculated frequencies of the trapped vertical p-mode oscillations for various parameters. We 
arrange the results from the viewpoint how the frequency change of the = oscillation and 
that of the rij- = 1 oscillation are correlated, when c\/c^ is changed, under other parameters 
being fixed. The correlation is examined to two different oscillation modes, i.e., n = l and n = 2. 
Figure 1 is for n = 1, while figure 2 is for n = 2. 

The correlation depends on M and a^. In figure 1, two cases of M = IAMq and M = 
1.8Mq are shown. The cases of M = 1.4Mq are shown by thick curves for two different values 
of spin parameter, i.e., a* = and a^, = 0.2. The thin curves are for M = I.SMq and two cases 
of a^, = and a^, = 0.2 are shown. The dots attached on each curves show the value of cj^/c^. 
When c^/Cg is small, both frequencies of the n,- = and = 1 oscillations are high and the 
point for the c^/Cg is on the upper- right corner of the correlation curve. As c^/Cg increases, 
the point move towards lower-left corner along the curve, as shown in the figure. It is noted 
that the correlation curves shift downward as M increases, while upward as a* increases. 

In figure 2, the frequency- frequency correlation is shown for oscillations with n = 2, as 
mentioned before. The set of (M, a^) adopted is, however, different from that in figure 1. That 
is, the thick curves are for M = IAMq as in figure 1, but the values of a* adopted are a* = 0, 
0.4, and 0.6. The thin curves are for M = 1.2Mq, and the value of a* adopted are = and 
0.4. 

One of major differences between oscillations of n = 1 and n = 2 is that if the value of 
c^/Cg is the same, the oscillations of n = 2 (figure 2) have lower frequencies compared with 
those of n = 1 (figure 1). Dependences of correlation curves on M and a* are qualitatively the 
same in both cases of = 1 (figure 1) and n = 2 (figure 2). 
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Fig. 1. The frequency correlation between tlie ~ and = 1 oscillations, when c\/c^ is changed. 
This figure is drawn for oscillations with n = l. The correlation depends on M and a*. The thick curves 
are for AI = IAMq and two cases of a* = and a* = 0.2 are shown. The thin curves are for M = I.SAIq 
and two cases of a* = and a* = 0.2 are shown. The dots attached on curves are the value of c^/cg . 




Frequency of nr= 1 Oscillations (Hz) 



Fig. 2. The frequency correlation between the = and = 1 oscillations, when c^/Cg is changed. 
This figure is drawn for oscillations with n = 2. The thick curves arc for M — IAMq and three cases of 
a* = 0, a* = 0.4, and a* = 0.6 are shown. The thin curves are for M = 1.2Mq and two cases of a* = and 
a* = 0.4 are shown. The dots attached on curves are the values of c\/cl as in figure 1. 
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5. Comparison with Observed Frequency Correlation 

The next problem is to compare the frequency-frequency correlation shown in figures 
1 and 2 with the observed correlation between the lower and upper kHz QPOs. The latter 
has been plotted on a frequency-frequency diagram for some typical neutron-star X-ray sources 
(e.g., Abramowicz 2005; Belloni et al. 2007). Here we adopt the figure by Abramowicz (2005) 
and superpose it on figures 1 and 2, which are shown respectively as figures 3 and 4. In order 
to avoid complexity, the values of c\/cl as well as M and a* are not shown in the superposed 
figures. 

Comparison of figures 3 and 4 shows that the correlation in the case oi n = l (figure 3) 
better fits to the observed correlation of twin kHz QPOs, compared with that of n = 2 (figure 
4). For example, in the case of n = 1, the observed correlation points of Sco X-1 (green points) 
are well on the correlation curve calculated by M = 1.4M0 and a^. = throughout the whole 
observational points. In this case of M and a*, the range of variation of c^/Cg necessary to 
describe the observations of Sco X-1 is roughly c\/c^ = 2 ~ 6. The set of (M, a*) that well 
describe observations are, however, not unique. In the case of Sco X-1, for example, we can 
describe observational data by adopting a higher mass than I.4M0 and a non-zero a*. That 
is, if we adopt M = I.SMq, for example, the curve of a* ~ 0.25 fits well the observational data 
(figure 3). 

One of reasons why the n=l oscillations will be more realistic than the n = 2 ones is that 
in the latter case a rather high spin is required to describe observations if we take M > I.4M0 
(figure 4). In other words, M < I.4M0 is required if the spin is taken to be moderately small, 
say a* < 0.2 (see figure 4). 

If we take the picture that the observed twin kHz QPOs are the n = 1 oscillations, a 
natural question raised is whether the n = 2 oscillations are observed or not. One of possibilities 
is that the n = 2 oscillations are not observed since their amplitude will be small. The n = 2 
oscillations have one more node in the vertical direction compared with the n = l oscillations. 
Hence, the saturated amplitude of the n = 2 oscillations may be smaller than that of the n = 1 
Another possibility is that the QPOs observed in a low frequency region (say, 150 
Hz ~ 300 Hz) are the n = 2 oscillations. This is because if we want to describe low frequency 
QPOs by the n = 1 oscillations, too strong magnetic fields are required (see figures 1 and 3). 

6. Summary and Discussion 

In this paper we take the standpoint that the kHz QPOs observed in neutron-star X-ray 
binaries are disk oscillations trapped in the innermost region of relativistic disks. There are 

^ The saturated amplitude will be determined by the balance between viscous damping of oscillations and 
stochastic excitation by turbulence. The viscous damping will become strong as n increases, since the vertical 
wavelength becomes short with increase of n. 
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Fig. 3. Diagram comparing the calculated frequency- frequency correlation with the observed frequency 
correlations in some typical neutron-star X-ray sources. The oscillations with n = 1 arc adopted. This 
figure is a superposition of figure 1 and the observational data plotted by Abramowicz (2005). The straight 
line labelled by 3:2 is the line on which the frequency ratio of twin QPOs is 3:2. The four sources on the 
line are black- hole candidates and are outside of our present concern. 




Fig. 4. The same as figure 3 except that n = 2 oscillations are considered here. That is, this is a 
superposition of figure 2 and the observational data plotted by Abramowicz (2005). 
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many kinds of trapped oscillations in the inner region of disks, if the disks have inner boundary. 
One of important checks examining which trapped oscillations can become a candidate of QPOs 
is to study whether they can describe the frequency correlation of observed twin QPOs. This 
problem, however, seems not to be examined quantitatively yet. In this paper, we examined 
this problem to the vertical p-mode oscillations. The results show that in the conventional 
values of parameters of mass and spin, the vertical p-mode oscillations seem to well describe 
the observed correlation of kHz QPOs. 

The magnetic fields required in this model, however, are rather strong, i.e., /5 defined by 
/3 =Pg/Pmag = / c\ is less than unity, where Pg and Pma.g are, respectively, the gas and magnetic 
pressures. It should be noted that the presence of steady low-/3 (i.e., (3 <1 or c\/c^ > 2) disks 
has been numerically demonstrated by Machida et al. (2006) and Oda et al. (2007, 2009, 2010). 
One of reasons why magnetically supported disks are possible is that the escape of magnetic 
fields from the disks by the Parker instability is suppressed by strong magnetic tension (Shibata 
et al. 1990). A low-/3 disk is a bridge between ADAF and slim disk, describing a bright hard 
state (optically thin) and a high/soft state (optically thick) (Oda et al. 2009, 2010). In their 
numerical models, such a high value of c\/c^ as cj^/c^ = 100 has been realized in extreme cases. 
As far as the geometrically thin approximation is concerned, such a high value of c\/c^ is 
acceptable. For example, for c\/c^ = 50, equation ([3]) gives H/r ~ Scg/rfi. In the standard 
disks, the value of Cs/rQ is about 1/30 in the inner region of disks (e.g., see figure 3.4 by Kato 
et al. 2008). In low-/3 disks, however, the temperature is higher than that of the standard disks 
by the presence of dissipation of strong magnetic fields, and is between the temperatures of 
ADAFs and of the standard disks. If we take a temperature 10 times higher than that of the 
standard disks, we have Cs/rQ ~ 1/10 and H/r ~ 0.5 for c^/Cg = 50. 

In this paper, we assume that the disk is isothermal in the vertical direction. If the 
vertical structure of the disks is not isothermal but polytropic, the polytropic index is an ad- 
ditional parameter contributing to frequencies of trapped oscillations (Kato 2010). That is, 
the frequencies of trapped oscillations depend not only on the strength of toroidal fields, but 
also on the polytropic index. The frequencies of trapped oscillations also depend on boundary 
conditions (see Kato 2011a). This suggests the possibility that we can get more detailed infor- 
mation of the innermost region of disks by detailed adjustment of calculated correlation curves 
to observational ones for individual sources, if the present model of QPOs is correct. 

It should be noted here that the trapped region of the vertical p-mode oscillations is 
rather narrow in the case of n = 1 and n = 2, i.e., the outer edge of the trapped region is 
around 4rg when a* = (see figure 1 by Kato 2011a). A problem is whether the oscillations 
in such a narrow region can produce observed luminosity variations. One of another problems 
related to the present model is that the magnetic fields required to describe the low frequency 
QPOs (100 Hz ~ 300Hz) is rather high (e.g., c\/c^ ~ 10^), if the QPOs are regarded as the 
n = 1 oscillations. This may suggest that the sources with such low frequency QPOs have really 
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strong magnetic fields or these QPOs are oscillations of n = 2. 

Whether trapped vertical p-mode oscillations are really excited in disks is a problem 
remained to be examined. We suppose that they will be excited by stochastic processes of 
turbulence in disks, as solar and stellar non-radial oscillations are now known to be excited by 
the processes in the convection zone. Turbulence has two effects on oscillations. One is their 
excitation by stochastic processes and the other is their damping by usual turbulent viscous 
processes. In the case of disks, the turbulence is MHD turbulence resulting from the magneto- 
rotational instability (MRI), and is much stronger than that in stellar convection zone. Hence, 
the turbulence will excit and maintain trapped oscillations at a high level against dissipation. 

Finally, high-frequency twin QPOs observed in black-hole X-ray sources are briefly men- 
tioned. These QPOs have robust frequencies with their frequency ratio 3:2. Hence, the present 
model cannot describe these QPOs, unless in the black-hole cases we have particular processes 
which keep the strength of magnetic fields unchanged with time by balance among creation (by 
differential rotation and MRI), escape (by jet and accretion), and dissipation (reconnection) of 
magnetic fields. In addition, the frequency ratio must be kept to 3:2. Considering these situa- 
tions, we think that some additional mechanisms will be necessary to describe the twin QPOs 
of black-hole X-ray binaries in the framework of the present model of trapped oscillations. 
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